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CNJ ' A pulsar model is proposed which involves the entire magnetosphere in the produc- 

tion of the observed coherent radio emission. The observationally-inferred regularity of 
^j' 1 peaks in the pulsar profiles of 'slow' pulsars (Rankin 1990,1993a), is shown to suggest 

that inner and outer cones of emission near the polar cap interact with and 'mirror' 
two rings in the outer magnetosphere: one where the null line intersects the light- 
cylinder, and another where it intersects the boundary of the corotating dead zone. 
The observed dependency of conal type on period is shown to follow naturally from 
CNj ■ the assumption that cones only form when the mirror intersection points lie between 

!> ' two fixed heights from the surface, suggesting that a feedback system exists between 

the surface and the mirror points, accomplished by a flow of charges of opposite sign 
in either direction. In their flow to and from the mirror points, the particles execute an 
azimuthal drift around the magnetic pole, thereby creating a ring of discrete 'emission 
nodes' close to the surface. Motion of the nodes is observed as subpulse 'drift', which 
■ is interpreted here as a small residual component of the real particle drift. The nodes 

\ can move in either direction or even remain stationary, and can differ in the inner and 

outer cones. A precise fit is found for the drifting subpulses of PSR0943+10. Azimuthal 
interactions between different regions of the magnetosphere depend on the angle be- 
tween the magnetic and rotation axes and influence the conal type, as observed. The 
Q ■ model sees 'slow' pulsars as being at the end of an evolutionary development where 

| the outer gap region no longer produces pair cascades, but is still the intermittent 

c/3 , source of low-energy pairs in a magnetosphere-wide feedback system. 

a ■ 

• • ■ 1 INTRODUCTION tional grounds too it is not easy to reconcile the model 
m ^ ' with the complexity of subpulse behaviour, such as star 
>v> , In 1975 Ruderman & Sutherland (henceforth RS) proposed tionary or counter-drifting subpulses (e.g. Biggs et al. 1985, 
» I ■ a model for pulsar emission which postulated a gap region Nowakowski 1991), mode changing and nulling. Evidence 
d located immediately above the neutron star surface and, as that the emission profiles generally take the form of nested 
a result of the inability of the electric field to remove ions cones (Rankin 1983, Lyne & Manchester 1988, Rankin 
from the surface, subject to an intense potential difference 1990(RI), 1993a(RII), 1993b) has presented a further chal- 
of around 10 14 eV. In this gap pair- production in the intense lenge to the model, and only by appealing to surface multi- 
magnetic field generated electrons which heated the sur- poles (e.g. Gil et al, 2002a, b) does it seem possible to confine 
face and created localised discharging regions. These 'sparks' the 'sparks' either radially or azimuthally. More recently, the 
drifted around the magnetic pole and ejected energetic par- precisely drifting emission columns of PSR0943+10, anal- 
ticles into the magnetosphere where a bunching mechanism, ysed in great detail by Deshpande & Rankin (1999, 2001 
first proposed by Sturrock (1971), created coherent radia- (henceforth DR)), convincingly confirm the picture of cir- 
tion following a two-stream instability. The model has been culating plasma columns. However the observed circulation 
highly influential, since it provides a quantitative framework rate is only reconcilable with the RS model by adopting 
within which theorists and observers can work, and one of its a potential difference across the cap which is significantly 
underlying tenets, that the sources of pulsar radio emission lower than that predicted by the RS model, and - if gener- 
are plasma columns circulating just above the polar cap, is ated by surface magnetic features - the columns would seem 
now widely accepted. to require multipoles with a high degree of regularity (Asseo 

v , , , , , & Khechinasvili 2001, Gil et al 2002a,b). 

Nevertheless, over the years the polar gap model has ' ' ' 

been questioned on theoretical grounds, largely through dif- In response to these difficulties, some authors have anal- 

ficulties with the neutron star surface binding energy (e.g. ysed magnetospheres based on free electron flows from the 

Jones 1985, 1986, Abraham & Shapiro, 1991, Neuhauser polar cap and gradual acceleration into the upper magne- 

et al, 1987) and with the high plasma densities required tosphere ('inner accelerator' models) (e.g. Arons & Scharle- 

in the emisson regions (Lesch et al, 1998). On observa- mann 1979, Mestel & Shibata 1994, Hibschmann & Arons 
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2001a,b (henceforth HAa,b), Harding & Muslimov 2002 
(HM), Harding et al 2002 (HMZ), Harding & Muslimov 
2003). Such models are often part of wider attempts at cre- 
ating a self-consistent global magnetosphere (Mestel et al 
1985, Shibata 1995, Mestel 1999 and references therein). Al- 
though arguably more consistent with the physics of the 
neutron star surface - and with the overall current balance 
and torque transfer requirements - these models lack the 
predictive power of the RS model when faced with detailed 
radio observations. Possibly the belief prevails that radio 
emission, undoubtedly originating just a few tens of stellar 
radio above the surface and energetically weak compared 
to the spin-down energy loss, has little to say about global 
conditions. Furthermore the sheer complexity of highly time- 
dependent subpulse phenomena acts as a great deterrent to 
theoreticians seeking to define a steady state condition. 

In trying to bring the large-scale magnetosphere models 
closer to observational testing, and in a deliberate attempt 
to explore the link between the the observed polar cap ac- 
tivity and the outer magnetosphere, we have taken a fresh 
look at the radio profiles analysed in RI and RII. Pulsar in- 
tegrated profiles are suitable to this purpose, since they are 
one feature of pulsar observations which displays great sta- 
bility: a pulsar's profile is its invariant signature. Rankin's 
conal classification of the profile forms is here interpreted es- 
sentially from a geometric standpoint without prejudice for 
or against any particular formation process. In essence, we 
seek to use the profiles as diagnostics of the magnetosphere's 
structure. 

This paper argues that if regular conal structures exist 
in all or most pulsar profiles - as claimed in RI and RII - 
then explanations in terms of complex polar magnetic field 
topology are difficult to support, since they would require 
similar magnetic topologies from pulsar to pulsar. However 
it is suggested here that conal emission may be the natural 
geometric result of interactions between the polar cap and 
outer regions of the outer magnetosphere in the purely dipo- 
lar environment of 'slow' older pulsars, provided that the 
once prolific outergap production of pairs of a pulsar in its 
early life is now weak and intermittent. This proposed evolu- 
tionary link between the younger faster pulsars (such as the 
Crab pulsar) , whose principle emission is in the gamma-ray 
band, and the older weaker radio pulsars is seen as a major 
and novel feature of the model here. 

First, in section (2), it is argued that the fixed ratio of 
the cone radii is consistent with the assumption that emis- 
sion occurs preferentially on two critical sets of field-lines, 
one bounding the corotating dead zone (defined by the last 
field- line to close within the light-cylinder), the other pass- 
ing through the intersection of the null line (defined as the 
boundary between regions of net charge density of oppos- 
ing sign) and the light cylinder. In section (3) it is shown 
that the observed period dependence of conal types can be 
simply explained if emission is possible only if a 'mirror' 
point (i.e. an intersection of the null line with the critical 
field-lines) lies between two fixed altitudes. In section (4) it 
is shown how the inevitable drift of the inflowing and out- 
flowing particles about the magnetic axis causes 'emission 
nodes' to form above the surface, and, in section (5), that 
this leads to their precession around the magnetic pole, gen- 
erating the well-known phenomenon of 'drifting subpulses' 
and hence the double cone structure. In section (6) the drift 



model is applied to the pulsar PSR0943+10 and precise fits 
to both the B-mode and the Q-mode are found. In section 
(7) the global nature of the pulsar phenomenon is stressed: 
azimuthal interactions between critical regions of the outer 
magnetosphere are shown to explain the observed depen- 
dency of conal formation on the angle of inclination. Finally, 
in section (8), the physical requirements of the underlying 
feedback system are discussed in the light of current physical 
ideas. 

The sketched model which emerges from this analysis 
contains many elements of the RS model, but on a scale 
which involves the entire magnetosphere: particles (pre- 
sumed here to be electrons, although the system is sign- 
reversible) are accelerated from low Lorentz factors close to 
the polar cap to achieve high 7 near the outer gap, which 
stretches from the corotating dead zone of the magneto- 
sphere to the light cylinder. Pair creation occurs in these re- 
gions, although certainly not in the profusion of the pair cas- 
cades supposed in Crab-like pulsars (Cheng et al, 1986, Ro- 
mani & Yadiaroglu, 1995). Rather the production is likely to 
be intermittent, weak (ie of low multiplicity with low Lorentz 
factors) and azimuthally-dependent, with most of the parti- 
cles produced forming a wind beyond the light cylinder, but 
with a small but essential fraction of the positrons returning 
to the surface, accelerated and funnelled by the increasingly 
tight bundle of magnetic field-lines above the poles (Michel 
1992). 

Unlike the case of fast pulsars (Cheng et al, 1986), the 
downward flow of particles is inadequate to screen the polar 
potential. Shortly before reaching the surface, the positrons 
emit sufficiently energetic radiation to cause a bunched 
avalanche of pairs which bombard the surface with coherent 
radiation. This is then reflected and contributes to the core 
component of the pulsar profile. Residual electrons formed 
by this process (maybe augmented by electrons emitted from 
the heated spot on the surface) are then accelerated back 
to the outer gap. This process is then repeated, and is re- 
inforced if the combined and equal drifts of the electrons 
and positrons around the pole maintain the hot spots at 
locations which are either fixed in azimuth or drift slowly 
backwards or forwards. The outflowing electrons are some- 
how stratified by the inflowing layers of pairs, and radiate 
curvature radiation parallel to the field lines and at about 
200 km (for 1GHz) above the surface. The holistic nature of 
the model means that the radio emission can simply be seen 
as an image of, and as driven by, the activities of the outer 
gap. 

The emission, particle flow and pair-creation processes 
are not worked out in detail here, and in many cases have 
been cannibalised from existing models (especially Mestel et 
al (1985), Michel (1992), Shibata (1994), HAa,b, HM, HMZ, 
Hirotani & Shibata 1999, 2001 (HSa,b)). Some features lack, 
as yet, a proper theoretical investigation: for example, the 
mechanism and level of weak low-energy pair production 
in outergaps has never been explored, since hitherto it has 
been assumed that outergaps only play a significant role in 
the emission of fast, young pulsars. On the other hand, the 
recent insight (see HAa,b, HM, HMZ, above) that in older 
pulsars, whose pair production at the polar cap relies on In- 
verse Compton Scattering, there will be a residual potential 
to accelerate particles up towards the outergap (and by im- 
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plication could accelerate returning particles of an opposite 
charge) gives support to a fundamental aspect the model. 

However we strongly stress througout the paper that 
the direction and purpose here is to establish the broad 
characteristics of a workable model, based on observational 
rather than theoretical grounds. In doing so we clarify not 
only what features such a model needs, but also what is not 
needed: there is no polar gap, no pair-creation in the outflow 
before the outer gap is reached, no self-stratification of the 
outflow to generate the coherence of radio emission, and no 
multipoles. 



2 INTEGRATED PROFILES 

From the early days of pulsar research it has been known 
that despite the complexity of individual subpulse behaviour 
the integrated profiles of pulsars remain remarkably stable 
(Taylor & Huguenin 1971). The profiles themselves adopt 
many forms, and despite many years of patient research 
(Rankin 1983, Lyne & Manchester 1988, RI, RII, Rankin 
1993b, Gil & Krawczyk 1996) no consistent morphology has 
yet found universal acceptance. Although alternative inter- 
pretations can be argued (e.g. the 'patchy' models of Lyne 
& Manchester 1988, Han & Manchester 2001, Smith 2003), 
we will here adopt the analyses of RI, RII, Gil et al (1993) 
and Kramer et al (1994), which suggest a picture which is 
strikingly simple: at any given frequency, any pulsar profile 
can be represented as a cut across a notional emission en- 
velope made up of a central core plus two concentric cones. 
The core is centred on the magnetic axis above the pulsar 
magnetic pole, and the cones surround it one or two hundred 
kilometres above the polar cap. Sometimes one or other of 
the cones is missing. As the frequency decreases the cones 
widen according to a dipole geometry (radius-to-frequency 
mapping), so that the conical structure is 'tied' to the field- 
lines. We propose to explore the double cone model in order 
to see what geometric - and hence physical - consequences 
flow from it. 

In Rankin's work (RI, RII) much depends on the ac- 
curacy of the inferred values for a, the angle between the 
pulsar's rotation and magnetic axes. These are calibrated 
using the measured angular width of the core component in 
perpendicularly rotating pulsars. This method implicitly as- 
sumes that the core component always spans the same field- 
lines, that at all alignments the core component appears to 
be radially emitted from the surface in the region of open 
field-lines between the boundaries of the closed 'dead' coro- 
tating zone. These assumptions have the obvious weakness, 
conceded by Rankin in RI, that no physical model has hith- 
erto been suggested with these properties. Yet a possible 
explanation for this (eg Michel 1992) emerges later from the 
analysis here in terms of reflected emission from downfalling 
particles. However, we stress from the start that our purpose 
is not to impose such a model, but merely to point out that 
the geometry will support it. 

The aim of the detailed analysis in RI and RII was to 
disentangle the geometric coincidence created by our line 
of sight from the underlying intrinsic geometry. Using the 
established values of a to infer the geometry of the inner 
and outer cones, it was found that, when present, the cones 
always subtend the same angular radii to the magnetic axis, 
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Figure 1. The geometry of the 'mirror' points in a plane pro- 
jected onto the rotation axis. The null line intersects the dead 
zone and the light cylinder at r* dz and r ( * c respectively. 



dependent only on the period P and frequency. In a study 
of some 150 pulsars Rankin (RII) concluded that the outer 
and inner radii appropriate for 1GHz are 

5 75 o p -0.50 



Pouter 

and 



pinri 



= 4.33°P~ 



(1) 



(2) 



The symmetry of the conal positions about the centre of 
the profile is all the more surprising since the intensities of 
the peaks are often highly asymmetric (e.g. PSR1133+17 
(Nowakowski 1996)). Some pulsars possess both cones (type 
M), others merely one (type T), and those with one cone 
may take either the inner or the outer radius. Type St, a 
young population generally exhibiting a single core compo- 
nent, sometimes reveal a conal structure at high frequencies 
(RII). It is remarkable that, whatever the configuration, in- 
termediate values between (1) and (2) are rarely found. Fur- 
thermore, most pulsars seem to have no more than two cones 
(although four or five relatively fast pulsars may have cones 
displaced onto wider cones, but in the same ratio (Mitra & 
Deshpande, 1999, see also Section 5.4)). 

The double-cone results were apparently confirmed 
(with smaller samples) by Gil et al (1993) at 1.4GHz and 
10.55GHz, and Kramer et al (1994) at 1.4GHz, 4.75GHz 
and 10.55GHz, yielding the expected smaller opening angles 
consistent with the radius-to-frequency mapping. None of 
these studies included millisecond pulsars, whose properties 
will also not be considered here, in a deliberate attempt to 
compare only pulsars with similar magnetic field strengths. 

In interpreting her results, Rankin (RII) notes that the 
inverse root dependence on period supports a dipole geome- 
try: the radius of the polar cap, defined by the last field-line 
to close within the light-cylinder, scales with P~ 0,5 , and she 
therefore suggests that the two cones are both formed tan- 
gential to this same field-line at differing heights, with the 
emission of the inner cone lying below that of the outer cone. 
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Following the logic of the dipole geometry, the outer cone at 
1 GHz, whose angular radius is given in (1), is then formed 
at a period-independent height of 220km above the surface, 
and the inner cone at 110km. 

This interpretation is difficult to reconcile with the 
frequency-to-radius mapping concept, whereby greater 
heights correspond to lower frequency emission, and theoret- 
ical efforts have been made to resolve this. Petrova (2000), 
Petrova & Lyubarskii (2000) and Qiao et al (2000) have sug- 
gested bimodal propagation models for the inner magneto- 
sphere. However these models have been devised to explain 
the conal profiles, and it is as yet hard to see how they 
operate at the more fundamental subpulse level without im- 
plying that the inner and outer cones have identical subpulse 
behaviour (something which is rarely, if ever, observed). 

In this paper we adopt an alternative interpretation: 
that the inner cone is also formed at the same height as the 
outer cone, but on a field-line closer to the magnetic axis. 
This suggestion is not new, and can be made consistent with 
the RS model by including surface multipole components of 
a specific form (Gil et al 1993, Gil & Sendyk 2000, Asseo & 
Khechinasvili 2002). But it would also enable the observa- 
tional conclusions of Kijak & Gil (1997) and Kijak (2001) 
that emission heights decrease with decreasing period to be 
reinterpreted as emission cones moving closer to the mag- 
netic axis. 

Clearly identifying the true fieldlines on which emission 
occurs is a formidable observational task (Mitra & Rankin 
2001, Kijak & Gil 2003). But in the present context, where 
we are pursuing the hypothesis of a purely dipolar field, any 
suggestion that different cones are located on different field- 
lines is highly radical, since it would appear to imply that 
special field-lines are somehow chosen near the surface by the 
upward-moving particles - and the question arises how the 
particles can 'know' which field-lines to select! The dilemma 
can only be resolved by abandoning the concept, long-held 
by theoreticians, that the structure of the magnetosphere is 
determined, indeed driven, by conditions near the polar cap. 
The step we take here is to reverse this logic, and argue that 
events in the outer magnetosphere are the true 'engine' of 
the pulsar system, and that the observed polar region events 
are passive reflectors of these. 

In their seminal work of 1969, Goldreich & Julian(GJ) 
divided the polar cap into a central region of outflowing 
(negative) particles and an outer annulus within which the 
current returned to the star. The outer boundary was de- 
fined, as in Rankin's work, by the last closed field line. The 
inner boundary was defined by the field-line which intersects 
the light cylinder at the same point as the so-called null sur- 
face. The null surface is to first approximation the locus of 
the points on which the magnetic field is perpendicular to 
the rotation axis (see Mestel 1999, p536), and is physically 
significant since it is the surface separating regions of net 
positive charge from those of net negative charge on which 
the net charge density must be zero in a steady system. In an 
aligned pulsar the central region has negative charge density 
and the annulus positive. This polarity is reversed in the the 
counter-aligned model of RS, where the polar cap does not 
extend beyond the inner boundary. 

The importance of the null surface was first realised by 
Holloway (1973), who pointed out that a charge-separated 
particle flow could never smoothly cross it, and that there- 
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Figure 2. The geometry of the 'mirror' points of a near- 
pcrpcndicular rotator in the plane of the observer's line of sight. 
The null line intersects the dead zone at as before, and the 
light cylinder at a slightly reduced r* c . 



fore an 'outer gap' of intense electric field would form, ex- 
tending from the corotating dead zone to the light cylinder 
surface. In energetic young pulsars this gap is thought to 
be the site of gamma radiation (Cheng et al., 1986). In an 
axisymmetric system the null surface forms a cone of half- 
angle arctan V2 ~ 54.7°, which intersects the light cylinder 
at an altitude of 



, 3 . 2 3 2 C P 

and the dead zone at 



2 cP 
S R ' C 3^ 



(3) 



(4) 



where Ri, 



Mc — if- is the light-cylinder radius. The field-lines 
linked to these extrema enter the neutron star surface at 
angles which are in the period-independent ratio of 



(^)~ =Q*=0.74 (5) 

(inner to outer). Through dipole scaling this ratio would 
then apply at any given height above the surface. 

Of course real pulsars are inclined rotators and in their 
outer magnetospheres are far from axisymmetric and far 
from dipolar. Yet our line of sight passes over the magnetic 
pole with a low 'impact' angle in a plane where the magnetic 
field close to the surface is dipolar and hence axisymmetric 
about the pole to a first approximation. In this plane the 
null lines (along which the magnetic field is perpendicular 
the rotation axis) remain unchanged from those shown in 
Fig 1 for a precise dipole, and will make the same intersec- 
tion with the dead zone at r* dz . However the light cylinder 
intersection will shift by an extent dependent on the angle 
of inclination. Fig 2 shows the most extreme case when the 
the magnetic and rotation axes are near-perpendicular, but 
still using a dipolar structure. The critical field-line for the 
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inner cone has moved slightly outwards and the ratio of 0.74 
has changed to 0.82. 

In the more realistic Deutsch (1959) vacuum solution 
the field-lines near the light cylinder are swept back and 
would give an asymmetric and curved null lines at large in- 
clinations. Furthermore the limiting field-lines defining the 
boundary of the dead zone will be perturbed, quite possibly 
in an asymmetric fashion. But we stress that for most pul- 
sars these effects may only be of second order, and the ratio 

f r* 1 

(-4 s -) 3 m ay not stray far from the 0.74 to 0.82 range, espe- 
cially since statistically the angles of inclination for pulsars 
cluster around 35° and are only exceptionally found at high 
inclinations (see Figs 2 and 4 of RI, and Section 7 of this 
paper). 

Taking the observed ratio of the inner to outer cone 
radii at 1GHz from (1) and (2) we obtain 



pouter 



0.75 



(6) 



The ratio is slightly higher at higher frequencies: 0.77 at 
1.4GHz (Gil et al 1993), 0.78 at 4.75GHz (Kramer et al 
1997). These ratios are consistent with the range indicated 
from the theoretical analysis above and suggest that we 
might identify the critical lines of GJ with the defining field- 
lines of the emission cones. This, in turn, points to an inti- 
mate relation between locations on the null surface - whose 
precise location may vary from pulsar to pulsar - and the 
pattern of emission above the polar cap. Somehow the outer 
magnetosphere and the polar regions mirror one another. 
This is an idea not present in polar gap models or space- 
charge-limited models for 'slow' pulsars (eg RS, HAa,b, HM, 
HMZ), yet it immediately suggests an evolutionary link to 
the outergap models for fast pulsars (HSa,b). Thus we may 
suppose that, although now weak, outergap pair production 
may still play a crucial role in the radio emission process. 



3 CONAL TYPE AND THE CRITICAL 
FIELD-LINES 

The principal result of RII is that the profile types (St, 
M, T) show a dependence on period, yet no dependence 
on surface-related parameters such as B s , the surface mag- 
netic field strength, or the acceleration parameter In 
general, inner cones are found to form at shorter periods 
while outer cones form at longer periods, with an overlap 
range where both cones are present. If we assume that the 
outer cones and inner cones are formed by the same physical 
process at either null-line intersection, then this result can 
be reproduced by a simple geometrical argument. All that 
is needed is the assertion that a cone only forms when its 
corresponding intersection point lies within a fixed, period- 
independent, range of altitudes. 

The simplest way to see this is to consider two pulsars, 
one of period Pi and a faster one of period P2 = (§) 5 Pi. 
Then from (3) and (4) the outer cone of the slower pulsar 
is formed on exactly the same field-line as the inner cone 
of the faster. And since the opening angle of the null sur- 
face is independent of period, the magnetic field strength 
and the distance of the intersection from the surface remain 
the same, although now the intersection is with the light- 
cylinder of the faster pulsar rather than the dead zone of 



the slower pulsar. Thus whether a particular field-line ap- 
pears as an inner or outer cone is determined solely by the 
pulsar period. This is completely independent of whatever 
physical criteria must be met in order to form a cone. 

With this insight it is possible to use the observed pe- 
riod dependency of the various profile types to determine 
the range of periods within which emission cones are cre- 
ated. Rewriting (5) , the distances, r ; * c and r* dz , of the inter- 
section points from the star (which are proportional to Ri c , 
and hence to the period) are in the ratio 



* V ry J 



,3 « 
'2' 



1.84 



(7) 



Let us then suppose that a cone only forms when the inter- 
section points fall within r* min and r^ax- For r* dz to be in 
this range, and hence for an outer cone to form, the pulsar's 
period must, from (4), satisfy 

c c 
Similarly, the period range for inner cone formation is, from 
(3), 



( 2 1 r*rnin 



< P <2(-)37r!^L 
3 c 



(9) 



The first range maps onto the second by the factor (7). As is 
shown in Fig 1, long-period pulsars with wide light-cylinders 
will be found to have only an outer cone because the light- 
cylinder intersection lies above r* max , and fast pulsars only 
an inner cone bcause the deadzone intersection is below 
r min> with pulsars with an intermediate range of periods 
having two cones. This is exactly as reported in RII. 

In RII Rankin finds that type St and type T with 
only an inner cone have a mean period around 0.5s, type 
M with two cones have a mean period of 0.86s, and types 
T and D (double) with only an outer cone have a mean 
period of about 1.25s. If we assume that cones are only 
formed when the null-line intersection point lies between 
r min = 20, OOOfem and r* max = 70, OOOfem, then the range 
for inner cone formation is 0.32s to 1.15s, and consequently 
0.6s to 2.15s for outer cones. This implies a double cone (M) 
range from 0.6s to 1.15s. These figures successfully repro- 
duce the observed mean values, although the scatter around 
the mean is fuzzy (see tables in Rankin 1993b), and other 
hitherto unconsidered factors may be involved. Later in this 
paper it is suggested that the degree of electric screening 
near the magnetic pole may be such a factor, resulting, for 
example, in core-dominant type St pulsars having a signifi- 
cantly lower range than the other categories. 

Note that it is the distance from the star which deter- 
mines the appearance or non-appearance of a cone, and not, 
for example, the ambient magnetic field strength. This would 
depend on B s , a factor which Rankin finds not to be statisti- 
cally significant. Assuming a surface field of B s = 10 12 G, the 
fixed range in radial distance translates into a poloidal mag- 
netic field range of around B m in = 100G to B ma x = 4000G, 
a factor of (3.5) ~ 40. Actual estimates of B s vary by a fac- 
tor comparable to this, and a criterion based only on field 
strength at the intersections could not reproduce the preci- 
sion of the observed range without showing dependence on 
B a . We are left with a criterion based on a tight range in 
intersection altitudes, or, equivalently, in star to null line 
communication times from 0.07s to 0.23s. 
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Figure 3. the vertical line (P1,P2,P3) show different possible 
positions for the light cylinder and their intersections with the 
null line (r* cl , '*;* c2 : r ;* c 3)' Each gives rise to a dead zone which 
intersects the null line at (r^ zl , r 2 z 2' r dz3^ respectively. Only 
r dz2' r fcl' r dz3' r lc2 ue w ithin the fixed range of heights indicated 
by r* ■ and r^ ax and therefore become mirror points 



We stress that as yet no particular model for the cre- 
ation of emission zones has been imposed. The whole pur- 
pose of this paper is to use the observations to define and 
constrain the emission model. What has been shown so far 
is that the observed dependence of profile morphology on 
period (and period alone) can be explained by a specifying 
a fixed and relatively narrow range of altitudes within which 
the 'mirror' site can fall. Thus the problem of morphology 
has been decoupled from the problem of devising an emis- 
sion model. Observations seem to be forcing us to consider 
pulsar radio emission as a global pan-magnetosphere phe- 
nomenon, involving intercommunication between the polar 
cap and the outer gap. The required features of an emission 
model based on this interpretation are discussed in the next 
section. 



4 THE FORMATION OF EMISSION NODES 

In the previous sections it is argued that the observations 
of radio profiles support the existence of a feedback system 
betwen the star and the mirror points. This implies a flow of 
charged particles in both directions, guided by the magnetic 
fieldlines. Unless the electric fields parallel to the magnetic 
fields are universally screened (we discuss in Section 7 recent 
work by Harding and coworkers (HM,HMZ), who indeed find 
that in older pulsars this potential cannot be fully screened) 
it inevitably follows that the streams of charged particles 
will undergo an azimuthal drift, irrespective of how they are 
produced at either end. The sense of the drift will be the 
same for particles of either sign, and independent of whether 
they are moving outwards or inwards. It is also independent 
of whether we are dealing with a pulsar with rotation axis 



and magnetic moment in the same sense (as in the 'inner 
accelerator' models) or with an antipulsar (as in RS), where 
the axes are in opposing directions. 

The picture we adopt is as follows: electrons are grad- 
ually accelerated outward on the critical mirror fieldlines 
from close to the surface to the the mirror points on the 
null line, where they initiate pair-production in the outer- 
gap. This production will not and cannot be prolific since 
the interacting photons will have low energy and will fail to 
initiate the massive pair cascades found in younger pulsars 
(Cheng et al 1986, HSa,b). Nevertheless sufficient low-energy 
positrons are then supposed to be available, either through 
pair-production or through interaction with another mag- 
netospheric region (as envisaged in the previous section) , or 
both, to cause a backflow of positive charge to the surface. 
The total current flow will be dominated by electrons, since 
in a steady state the net charge density must always be neg- 
ative and close to the charge-separated GJ value. The drift 
will cause the positrons to return to the surface on a different 
fieldline from which the electrons came (see Fig 4), and they 
will be accelerated towards a different azimuthal location on 
the polar cap. Close to the surface the positron energies will 
be sufficient to trigger a pair-production 'avalanche' (as de- 
scribed in Michel 1992), which will cause electrons, either 
from the surface or trailing the avalanche (or both) to be 
accelerated upward into the magnetosphere. These will, in 
turn, continue to drift on their passage to the mirror re- 
gions. If theory can demonstrate that no particles from the 
surface are needed or even possible in this scenario, the en- 
tire model becomes sign-reversible and can apply to pulsars 
and antipulsars alike. 

The central proposal here is that in a quasi-steady state 
the avalanches above the surface will be confined to discrete 
'emission' nodes, arranged so that the flow forms a continu- 
ous yet finite stream of particles linking all the nodes (Fig 4). 
It is presumed that by remaining separated they avoid mu- 
tual interference and maintain conditions for pair-creation 
and the observed coherent radiation. In this section we de- 
scribe and quantify the process by which these nodes are 
formed, and in the next section we consider how observa- 
tions of subpulse drift can be interpreted as a movement of 
these nodes and be used as diagnostics of the conditions in 
a pulsar's magnetosphere. 

In the previous section it is argued that observations 
suggest that a high level of interconnectivity between the 
poles is a prerequisite for pulsar emission. Thus we may ex- 
pect the particles created at the mirror points will feed both 
poles and ensure that particle drift at one pole is coordi- 
nated with that at the other. This supports a view that the 
mirror points are far from being the passive reflectors of 
events close to the surface, but rather drive the entire drift- 
ing phenomenon by means of null line interactions between 
the cones, and possibly through inter-pole links. Here we 
analyse an axisymmetric system at a single pole, although 
the ultimate intention is to apply the resulting estimates in 
the context of an inclined interactive magnetosphere. 

The particle drift rate is estimated by an argument anal- 
ogous to that of RS, but applied to an entire axisymmet- 
ric magnetosphere rather than confined to a polar cap. In- 
evitably the estimates will be less precise, but nonetheless 
useful. We consider first a configuration with an outer cone 
only (i.e. with no inner cone screening). If there is a poten- 
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tial drop of AV along the central magnetic axis from the 
surface to a height at the level of the dead zone intersection 
(see Fig 1), then the typical potential difference available 
over the cylindrical radius of R dz between the axis and crit- 
ical fieldlines close to the deadzone is — AV, since ^there is a 
zero potential drop along field-lines bounding the dead zone 
and across the surface of the polar cap. Then taking the 
typical scale length of the electric field in the region of the 
dead zone as Rdz we obtain j^- as a measure of the mean 
non-corotational electric field perpendicular to the magnetic 
field B dz at a cylindrical radius of Rdz- Using this we can 
estimate the mean magnetospheric drift rate relative to the 
corotating frame, and in the opposite sense to the star's ro- 
tation, as 



1 | AEXB 

R W C 
cAV 2tt 



AV 



(10) 



BdzRdz P (AVmax)out 

where AE is the variation in the electric field from the coro- 
tational value, and 

2$ OU i 27rB dz R 2 dz 



(11) 



cP cP 

is the maximum potential available in the vacuum case, and 
&out is the magnetic flux through the polar cap extending to 
its dead zone boundary (Sturrock 1971). Equation (10) gives 
zero net drift in the inertial frame when = (AV ma x)out, 
and corotation if AV = 0. Note that (flD) out is formally 
negative. The net particle drift rate in the inertial frame 
((ilo)out = 75- + (^d)ou*) is thus proportional to the level of 
screening within the magnetosphere, a quantity which may 
well be time-dependent. Henceforth, primed quantities will 
indicate that these are measured in the corotating frame. 

An equivalent estimate can be made for the inner cone 
drift, although this will be more approximate since the field- 
lines will be swept back with a significant toroidal compo- 
nent, and particles are likely to have sufficient energy to 
leave the field-lines as they cross the light cylinder (Mestel 
et al, 1985). Assuming the potential difference of AV along 
the central axis now extends to a height at the level of the 
light cylinder/ null line intersection yields a potential differ- 
ence — AV between the magnetic axis and the light cylinder 
at Ri c . This effectively assumes a zero potential drop along 
the null line between r* dz and r ; * c . In a more realistic model 
(i.e. that of Mestel et al 1985 or Shibata 1990) this potential 
drop is not zero and drives the closure of the current loop 
outside the light cylinder. Using now Ri c as the scale length, 
this will produce the parallel formula to (10), namely 

I ( il D)in |~ p p2 - "77 7X77 r - ^ 12 > 



BicRl 



with (AVmax)in calculated from the magnetic flux 
through the inner cone: 



2$ in nB lc Rf c 



(13) 



cP cP 

$ in is less than the full polar cap flux $ out by a factor of 
r -4f- , which from (5) gives 

(n D )' in = ^(n D )out = d)Hn D ) ut = i.84(n D )^ ( (14) 

r dz £ 



Thus in the simplified magnetospheric model used here the 
particle drift rate in the inner cone is nearly twice as fast as 
that in the outer cone, if either is present alone. 

Having established the parameters of the particle drift, 
we are now in a position to explore how many emission nodes 
are created. The inflow/outflow process in the outer cone 
will, in the corotating frame, circulate once around the cap 
It is conjectured that the circulation sets 



m a time -7- 

up regular cycle, whereby the flow 'visits' p' nodes within q' 
turns and then repeats itself continuously (Fig 4) . It follows 
that in each turn the flow makes descents to the cap, 
where 



P 



2tt 



3tt 2tt 



1 



1' I M out I 2r* 



2 P I (Cl D )' mt 

3lT (AVmax) out 



2 AV 

For the inner cone the equivalent calculation is 
p ^ 2tv c ,2a 2ir 1 



(15) 



,2,1 2tt 

(_)2 7r 

V 3 J p I (n D ) , 



47T (AVmax)out 

~9 AV 



(16) 



The relations between n and (Qd) in the first lines of (15) 
and (16) are purely geometric relations, showing how |t- de- 
pends on the drift rotation rate, the distance to the nodes 
from the mirror points and the period of the pulsar in the 
corotating frame. As will be seen in the next section, they 
enable us to set lower limits on £7 in both cones, since in nei- 

ther cone can (S7u) exceed the pulsar's own rotation speed 
and counter-rotate (Ruderman, 1976). 

The second lines of equations (15) and (16) express how 
£7 measures the extent to which the optimum potential has 
been screened: the lower AV falls the larger £7 will become. 

To compute |t- we need estimates of (AVmax)out and AV. 
(AVmax) out is obtained from (11): 

(AV ma x)out « 1.35 x ^10 13 eU (17) 

where B12 is the surface magnetic field in units of 10 12 G. 
For AV we can take a normalising value of, say, AV = 5 x 
10 12 , which corresponds to the potential needed to accelerate 
electrons (positrons) to approximately 7 = 9 x 10 6 at which 
radiated gamma-rays may trigger pair-creation. However it 
must be stressed that the intention of the model here is not 
to impose a value for AV or 7, but to find a way to measure 
them. Inserting (17) and the estimate for AV we obtain from 
(15) and (16) 

7 "p^^v - } (18) 

for the outer cone, and 

^A 5x10 ' 2 ) (19) 
q> P 2 1 AV ' 1 ' 

for the inner cone. 

We can now return to our result of Section (3): that 
emission cones tend only to form when their mirror points lie 
between two fixed, period-independent altitudes. With the 
model here, the upper limit may be interpreted as a maxi- 
mum distance of the mirror points above the emission nodes 
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Figure 4. A view of the particle drift looking down onto the 
magnetic pole in the corotating frame. The motion of particles 
over distances of more than 20,000km to and from the 'mirror' 
points on the null line creates emission nodes immediately above 
the polar cap. The nodes, and their corresponding mirror points, 
may be stationary, or drift in either direction. The pattern shown 
corresponds to the simplest mode (i.e. all the nodes are created 
in a single turn). 



at which the required coherence of the emission nodes can 
be maintained close to the surface - a high degree of angular 
precision must be required to confine the node to its partic- 
ular (possibly drifting) azimuth and prevent its 'smearing'. 
At the other extreme, a too-dense clustering of nodes may 
cause interference between the electron and positron flows of 
adjacent nodes, preventing a steady flow between the surface 
and the mirror points. 



5 DRIFTING SUBPULSES 
5.1 Introduction 

The emission of many pulsars exhibits a systematic modu- 
lation known as drifting subpulses, whereby a succession of 
subpulses, identified here as the emission of adjacent nodes, 
appear in the pulse window and each gradually 'drifts' in the 
same sense across the window to the position of his neigh- 
bour within a timespan of P3 rotation periods. The drift may 
be fast (« 2 periods) or slow (> 10 periods), and the pat- 
tern repeats itself, often for long stretches, before sometimes 
switching to a different 'mode' of emission. 

In interpreting observations of drifting subpulses, the 
model suggested here shares with RS the picture of emission 
columns circulating around the magnetic pole. However the 
underlying flow which produces these columns, illustrated 
in Fig 4, is radically different and has greater flexibility. A 
steady system of stable emission nodes is created and main- 
tained near the surface by a continuous, or quasi-continuous, 
flow of particles to and from the mirror points. Depending on 
the rate of particle drift, which in the corotating frame can 
have only one sense, the nodes may be stationary, or precess 



positively or negatively around the magnetic pole. An equal 
number of pair-producing regions on the outer gap, either at 
the dead zone boundary (rj z ) or at the light cylinder (r ; * c ), 
will precess in tandem (Fig 4). In theory a very wide range 
of emission systems can be generated in this manner, some 
highly complex and chaotic. However the 'mode' selected by 
the pulsar will reflect the value of . 

There is a fascinating everyday analogy to the process 
suggested here. If a camcorder linked to a television set is 
pointed at the screen in a completely darkened room, and 
if an initial instantaneous flash (such as the striking of a 
match) occurs between the screen and the camcorder, then 
the image of the flash appears on the screen and feeds back 
into the camcorder. The image on the screen will then con- 
stantly change, creating patterns which may be regular or 
chaotic, depending on the angle at which the camcorder is 
held about its axis. It is a classical demonstration of how 
complexity can arise even in a simple feedback system: com- 
plexity of outcome does not imply complex input. 

In the previous section we estimated the particle drift- 
rate in the corotating frame. This is the appropriate frame 
for the calculation of the driving electric field, but the ob- 
served number and drift of the nodes depend critically on 
the observer's frame of reference. For example, if the ob- 
server were rotating at a rate less than, but in the same 
sense as, the corotating drift-rate of the particles, their ap- 
parent drift relative to the observer would be smaller than 
in the corotation frame, and he would see a larger number 
of nodes for every rotation in his frame. These nodes would 
remain in the same sequence as in the corotating frame, but 
would begin to repeat before a single turn was completed. If 
the observer moved counter to the drift, he would see fewer 
nodes in a single turn, but on the second and later turns the 
nodes would appear to overlap, and the order of the nodes 
round the axis would therefore be changed. This is an im- 
portant point, since the pulsar observer in his inertial frame 
will always turn counter to the sense of the particle drift 
in the corotating frame, and may well observe a node pat- 
tern very different from that in the corotating frame (as will 
be seen in the case of PSR0943+10, discussed in the next 
section). 



5.2 Critical parameters 

We consider a distant observer's view of the drift, and relate 
the underlying particle drift in his inertial frame to the num- 
ber of nodes and to the observable quantity P3. We begin by 
adapting the estimates in (15) and (16) to the inertial frame 

drift {Q.D)out/in = + (^o) out /in> but still assuming that 
the pulsar is aligned. In this frame the observer sees p nodes 
created in q turns by successive inflows and outflows (where 
p and q are integers with no common factors) so that every 
q th node is sequentially linked as the particle flow circles 
the pole. Then 



p 

q 



2tt 



(fio)c 



2r* 



for the outer cone, and 



p 

q 



2tt 



{Q D )in I 2r* 



(20) 



(21) 
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for the inner cone. Note that | and ^7 in the notation of the 
previous section are equivalent parameters in the inertial 
and corotating frames respectively. 

Each of the particle path distances r dz and r* c defines a 
mean transit time (jdz or t; c in units of P, the pulsar rotation 
period) from the nodes to the mirror points and back in the 
outer and inner cone respectively, so that 



Tdz = 



Pc 



— = 0.212 

3tt 



for the outer cone, and for the inner cone 
2r,* 



,3,1 1 



0.395 



(22) 



(23) 



where the estimates are those for the aligned rotator as in 
(3) and (4). The inner cone value is almost double that of 
the outer cone, reflecting the greater distance the flow has 
to travel to the light cylinder intersection. 

In practice, q will be taken as low (< 10) and any more 
complex flow will be seen as a modulation of a basic (p,q) 
pattern. Since the particle drift in the inertial frame can- 
not exceed the rotation speed (Ruderman 1976), we demand 
that I (tt D ) out/in |< Hence from(20) and (21) 



P 1 37T 

- > ~ -7T =4.7 

q r dz 2 

for the outer cone, and 

5 > J- « (^4, = 2.6 
q Tic 3 

for the inner cone. 



(24) 



(25) 



5.3 Conditions for steady subpulse drift 

To turn the approximate results (20) and (21) into exact 
equations we need to incorporate into {ClD)out/in a term 
which takes account of the rate at which the nodes move 
around the magnetic axis. After a total time of P3 = pPzP 
(following the notation of RS for P3) the entire ring of emis- 
sion nodes will have presented itself to the observer. Thus 
in a smoothly-flowing system (20) can be rewritten as 



. 2n q 2n 2tt q , 1 1 . 

TdzP P 3 P pr dz qPs 



(26) 



In the equivalent expression for (fto)in, n c replaces Tdz- As 
long as j Pi |> 1, the first term in the brackets of (26) 
always dominates over the second. Hence the drift will be 
only a fraction of the true particle drift. In general, there 
is no preference for the sign of P3, and the observed drift 
may appear to be in the sense of rotation or the reverse. 
When both inner and outer cones are present, their nodes 
may conceivably have differing drift directions. Opposing 
drifts at differing longitudes of the pulse window have been 
observed in PSR0540+23 (Nowakowski, 1991). 

The picture implied by (26) is of a simple integrated sys- 
tem with a near-continuous even flow of particles streaming 
from each node up to the mirror points and back, succes- 
sively visiting all of the drifting nodes in the system cycle 
time (evaluated in (28) below) before returning to meet the 
node from which it set out. More realistically, the flow may 
be made up of 'packets', possibly separated by the node- 
mirror point travel time (« O.lPsecs), but sufficiently con- 
tinuous for us to 'see' each node at our sampling rate of P 



seconds. Should this occur the observed emission could be 
separately modulated by phase and by intensity. 

In practice, observers will be able to measure P3 
(though not necessarily its intrinsic sign), possibly also p 
(as in PSR0943+10), and will wish to infer q and T dz /i c - 
From the analysis so far they will have at hand conditions 
(24/25), and equation (26). But there is one further useful 
constraint which links the known and unknown quantities. 
This essentially geometric point arises because the magni- 
tude of the residual angle drift, represented by the last term 
in (26), must be small enough (i.e. | P3 j large enough) to 
ensure that the system has just p nodes and does not accel- 
erate/decelerate to a (p-1) or (p+1) system (depending on 
the sign of P3). In other words, for a given p, q and T dz /i c 
there is a critical | P3 | below which a steady system cannot 
form. 

In establishing this condition we will examine in greater 
detail, in equations (27) to (34), the geometric and physical 
features of a near-steady conal system which underlie the 
result (32). The following equations apply to either cone, 
but are illustrated using outer cone parameters. 

Consider particles being emitted simultaneously from 
all p nodes. The flow from each node will then call at its next 
successive (not necessarily adjacent) node after just r dz ~ 
0.212 periods, the fundamental unit of time for the outer 
cone subpulse system of the pulsar, which will in general 
depend on both P and a in an non-aligned system. In this 
time the node will appear to advance by 



2tt_ 



-Tdz 



3pP 3 



radians 



76.4° 
pP A 



(27) 



The node will drift smoothly if this angle is small compared 
with the angular extension of the node. The flow will then 
continue to the next node, and by the time the round trip 
is complete, all p nodes will have been visited in the system 
cycle time 

t c = PTdz « 0.212pPsecs (28) 
and the flow will have traversed a total angular distance of 

2n 2tt 4 

2nq + —t c = 2nq + —T dz « 2nq + — (29) 

P3 "3 O.T3 

This result can be equivalently obtained from (26) by mul- 
tiplying (Q, D )out by t c - 

The residual final term in (29) (which we will refer to 
as the cycle drift angle) is given by 



_ 2tt jt_ _ 76.4° 

c "P 3 T&Ri 3P3 ~ P 3 



(30) 



and represents the net angle through which the system has 
drifted in a single cycle. It is a useful observable parameter in 
determining the characteristics of the system, and it enables 
us to determine how many cycles, N c , are required to rotate 
the entire system through a full circle, namely 

Pa_.3.P 3 P 3 



2?r 

Oc T d 



0.212 



N c must be an integer if the system is to be closed and 
have a simple repeatability. Then any 'streaky' features in 
the underlying flow will share the P3 periodicity, the phase 
periodicity of the nodes, since 



P 3 = N c t c 



(32) 
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For the inner cone the net drift angle, O c , is significantly 
larger (±|^) and iV c « 

A further significance of the cycle drift angle is that it 
can determine, for a given P3, how close the p-node system 
is to going over into anp+1 or p — 1 system. As P3 reduces, 
the cycle drift angle grows (see (30)), and if the system is to 
keep its identity as a rotating p-node pattern, 6 C must not 
exceed the separation (-^ry) between successive nodes in the 
p + 1 or p — 1 system (depending on the sign of P 3 ). Using 
(31), a minimum P3 for a given p is therefore set by 

I Pa I^M^HCEi!) = o,12Mi) (33) 
q 3nq q 

where the positive sign applies if P3 is negative, or more 
simply, 

N c > 2^1 (34) 

Combined with the conditions (24)/ (25), (33) and its inner 
cone equivalent enable useful limits on p and q to be set. 
It is immediately clear that for simple systems with q=l, a 
high value for p will require a high P3 (i.e. a low drift-rate): 
for example, if p=20 in an outer cone geometry, P3 must 
exceed 4. This point is particularly relevant to PSR0943+10, 
discussed in the next section. 

Using his estimate for (p,q) and taking a plausible value 
for r dz , the observer could now evaluate the drift-rate of the 
particles, and hence infer the potential difference AV from 
(18). But in pulsars which exhibit great regularity in their 
subpulse drift it is possible to take one further step. Note 
that a particularly simple class of drift states can exist if the 
pulsar adopts a P3 such that the value of the square bracket 
in (26) is an integer or a simple fraction. This would result 
in a harmonic relation between the underlying particle drift 
and the star's rotation rate, a factor which would make the 
flow steady, sustainable and hence conducive to the forma- 
tion of emission nodes. It may even be a precondition for a 
pulsar to be observable. In an outer cone, — — is 4.7 in the 
aligned dipole case (from (22)), and the mean value of this 
term will not differ greatly at inclinations with low a. This 
suggests that in regularly drifting pulsars the value of the 
square brackets is 5, and that the harmonic relation to the 
rotation rate would be especially clear if p were a multiple 
of 5. Furthermore P3 has to be positive in nearly-aligned 
pulsars. From (26) and (31) the harmonic condition implies 
values for P3 and Tdz of 

P3 = l(N c +±) (35) 

T<u = i(l + 4r) ( 36 ) 
5 qN c 

Pulsars which exhibit highly regular drifting subpulses may 
be expected to have parameters which satisfy these equa- 
tions. 

5.4 Mode-changing 

Many well-known pulsars maintain a steady or near-steady 
drift behaviour over many periods (even over thousands of 
periods in the case of PSR0943+10 discussed in the next 
section) . These pulsars appear to have low values of the in- 
clination angle a (RI) , and hence the axisymmetric model of 



the previous section would seem a reasonable approximation 
of their conditions. Examples are PSR0031-09 (Vivekanand 
& Joshi 1997, Wright & Fowler 1981b), PSR1944+17 (De- 
ich et al. 1986), PSR2319+60 (Wright & Fowler 1981a), PSR 
1918+19 (Hankins & Wolszczan 1987), and, most recently, 
PSR0809+74 (Lyne & Ashworth 1983, van Leeuwen et al 
2002). 

However each of these pulsars also has at least one al- 
ternative - and equally stable - drift mode, so a number of 
steady solutions to equation (26) must exist in a single pul- 
sar. By inserting integral values of N c and q into (35) and 
(36) an infinite series of discrete possible stable harmonic 
states can be generated. Which of these is selected by a par- 
ticular pulsar - at a particular time - may be determined 
by allowed values of £Id and Tdz- These are likely, in turn, 
to depend on the current flow generated from the emission 
nodes close to the polar cap and on the precise height of the 
null line above them (which are interrelated parameters ac- 
cording to the recent outergap model of HSa,b). In the outer 
cone, particularly simple harmonic ratios (26) between fin 
and the rotation of the star (such as |,|,|,|, §, etc) may 
only be attainable if if p is a low multiple of 5 and not a 
multiple of q. A mode change to another harmonic state can 
occur providing N c , p and q satisfy the constraints (24) and 
(34). A change in N c will cause a change in P3, a change in p 
will cause a change in P2 , the observed frequency-dependent 
subpulse separation, and a change in ^ will change the har- 
monic ratio (26). 

Mode changes may (from (36)) also be accompanied 
by changes in Tdz, the particle travel time between the 
nodes and the mirror points and back, therefore requiring a 
shift in the emission region. Studies of the mode changes in 
PSR0031-07 by Vivekanand & Joshi (1997) and Wright & 
Fowler (1981b) suggest a progressive shrinking of the profile 
from mode to mode, which may support this interpretation. 
More recent work by van Leeuwen et al (2003) argues for 
a similar effect in PSR0809+74. However in all the pulsars 
identified above the drift repetition rate, P3, although stable 
for a while, shortens abruptly and successively through two 
or three discrete values before returning to the highest P3 in 
a quasi-cyclical manner. A possible model which emerges for 
this class of pulsars is one of progressively decreasing levels 
of screening and build-up of potential, followed by a sudden 
return to high screening, and a cycle is established. 

The quasi-stability of the general chaotic systems sug- 
gests that, once achieved, the pulsar's potential will remain 
in such a state, subject only to a gradual secular change, 
then return to either the original state or some 'closer' al- 
ternative mode. Key to the stability of the subpulse be- 
haviour during these transitions is the speed with which 
the nodes reconfigure themselves following a change: if the 
relaxation timescale is too long, the pulsar emission may 
never settle down to a steady pattern. The intermittent, 
but non-periodic, activity observed in the the core regions 
of many pulsars (Seiradakis et al. 2000), whatever its ori- 
gin, is evidence that most magnetospheres are indeed in a 
quasi-chaotic state. 
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6 A MODEL FOR PSR0943+10 
6.1 Background 

PSR0943+10, with a period of l.lsec and B 12 = 2 (see DR), 
has very typical pulsar parameters, but exhibits an extraor- 
dinarily precise alternating drifting pattern (in the domi- 
nant B mode) and occasionally switches to an apparently 
disordered Q mode (Suleymanova & Izvekova, 1984). It is 
of great interest here partly because the inferred angle of 
inclination between the magnetic and rotating axes is just 
11° — 15°, suggesting comparison with the aligned model 
used to estimate paramters in the previous section. But it 
is also of interest since Deshpande and Rankin (2000) (DR), 
in their recent detailed study of this pulsar, have provided 
the only example so far where N, the total number of the 
emission nodes, has been directly deduced from a precise 
measurement of the circulation rate. The exceptionally sta- 
ble subpulse pattern with P 3o 6 S = 1.867 is shown to have 
exactly N=20 drifting nodes, with the observed drift in the 
same sense as the presumed E x B particle drift. 

An interpretation in terms of the present model is some- 
what hampered by the fact that the line of sight trajectory of 
this pulsar is too oblique for observers to determine conclu- 
sively whether we are seeing the geometry of an inner or an 
outer cone (DR give two alternative good fits with a = 11.6° 
for an inner cone, and a — 15.4° for an outer cone). However 
the facts that Pg is positive and that the pulsar is almost 
aligned suggest (from the discussion in Section 5) that the 
emission is from an outer cone. Since N is known (and as- 
suming initially that we are seeing a single system, so that 
N=p), we can immediately estimate the particle drift-rate 
from (20) and (22) as 

(fiz>)«.t « ^| (37) 

It can be immediately seen (and is formally expressed in 
(24)) that q < 4 if (Qo)out is to be less than the pulsar 
rotation rate. 



6.2 The B mode 

The relatively small value of the observed P3 enables us to 
further constrain the possible configurations. Firstly, for ei- 
ther type of cone we can see that for p=20, q=l neither 
the outer cone condition (33) nor its inner cone equivalent 
is met, so, in the observer's frame, all 20 nodes cannot be 
linked in a single traverse round the polar cap. Systems with 
an even value of q can also be excluded since they would re- 
quire an odd number of nodes. Thus q=3 remains the only 
possibility if all nodes belong to the same system. However 
it is theoretically possible that we are observing 2 indepen- 
dent but interlocking systems of 10 nodes each, or 4 of 5 
nodes etc. These systems cannot satisfy the criterion for 
q=l, but are possibilities with q=3. Nonetheless, the sys- 
tem's observed clockwork regularity argues powerfully for a 
single integrated system, and therefore the most likely can- 
didate is one with p=20, q=3. 

Accepting this as our working hypothesis, what charac- 
teristics will the system possess? It will connect every third 
node at time interval (from (22)) of approximately 0.212P 
until all 20 nodes are visited, taking t c w 20 x 0.212 = 4.24 
periods for the cycle. Assuming only that P3 is in the region 



of 1.85, then from (30) the net angular drift of the system 
after the time t c will be approximately 



76.4° 
1.85 



= 41.3° 



(38) 



implying a net drift of about 2° per node visit. The nodes are 
18° apart, so node visits are separated by 3 x 18 + 2 = 56°. 

To make the model precise we need more exact values 
for both P3 and Tdz- Noting that the net drift angle 6 C is 
close to 40°, gives from (31) an estimate of N c as 8.73, close 
to the integer 9. Let us therefore suppose that in reality N c 
is exactly nine, giving the flow an exact repeatability. Then 
the further conditions (35/36), which would apply if the flow 
is harmonically coupled to the star's rotation predicts P3 and 
Tdz to be 



P3 = |(9 + J) = 1.86666... 
6 



Tdz 



= ¥ 1 + sh>' 



= 0.2074 



(39) 
(40) 



(39) is precisely the observed value of P3 and is the only value 
between 1.66 and 2.00 which can give harmonic coupling 
(these values corresponding to N c at 8 or 10). 

From the the value of Tdz in (40) we obtain an adjusted 
cycle time of t c = 20 x 0.2074P = 4.15P, exactly one ninth 
of the observed P3 = 37.35. With the deduced values for P3 
and Tdz it follows that the particle drift in the observer's 
(inertial) frame has an angular speed of 



2tt 3 1 
P X 20^0.2074 



+ 



3 x 1.867 J 



2tt 



x 0.7500 



(41) 



so the particle drift is exactly three quarters of the star's 
rotation rate. This suggests that the pulsar has adjusted 
both the path length of its flow and its potential field so 
that this coupling occurs. 

Thus we are left with a model of extraordinary har- 
mony: P3 is 9 times the basic timescale of 0.2074P, the par- 
ticle drift circulates round the magnetic axis 3 times in ex- 
actly 4 periods, requiring 20 x 0.2074P (= 4.15P) to com- 
plete the 20-node cycle, which then repeats itself exactly 9 
times before returning to its starting point. 

In the corotating frame the picture is even simpler. The 
angular speed of the the particles (fii>) out will now be ex- 
actly one quarter of the pulsar rotation rate, and in the 
opposite sense. The equivalent equation to (41), relating the 
particle drift to the p' and q' seen in the corotating frame 
will be 



2tt 



P 20 L 0. 2074 1 x (3 x 1.867) ' 
— x 0.2500 



(42) 



It is clear that p remains unchanged from p at 20 but now 
q' = 1, whereas q=3. Thus in this frame all the nodes are 
visited consecutively in a single turn, with the cycle still 
taking 4.15 periods, but an observer would see a much slower 
node drift P 3 = 3 x 1.867 = 5.600. This finally enables us 
to estimate the underlying potential difference required to 
support this drift. From (10) it follows that three quarters 
of the available electric potential is screened, thus leaving a 
potential (from (18) with B 12 = 2, P=l.l, j£ = p' = 20) 
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of AV = 5.3 x 10 12 eV, enough to accelerate particles to 
7^ 10 7 . 



6.3 The Q mode 

For the irregular Q-mode, which comes on abruptly and 
shatters the clockwork regularity of the B-mode, an inter- 
pretation is more difficult. But a number of features stand 
out. DR note that at the outset there is a roughly 4P inten- 
sity fluctuation present, which could be the residual of the 
underlying particle drift periodicity of the 20 node system. 
This periodicity is never observed in the smoothly-flowing 
B-mode, whose variations are those of phase rather than in- 
tensity. We would interpret this as the particle flow between 
the surface and the nodes becoming disordered, though pos- 
sibly even continuing to drift at a rate similar to that in 
the B-mode. The particles now flow in 'packets' rather than 
streams, and give rise to the 'streaky' emission which DR 
report. 

But why does the flow become disordered? The most 
prominent feature in the Q-mode power spectrum, spread 
across most of the pulse window, equivalent to a periodicity 
of P3 = 1.293, arises from occasional periodic bursts against 
a generally chaotic background. But with p=20 this drift no 
longer satisfies the condition (33) for any integral value of q. 
The possible suggestion is that the Q-mode has become dis- 
ordered because P3 has suddenly changed to a value which 
destabilises the 20-node pattern. Substituting the new P3 
into (30) and (31) yield a net drift angle, O c , of 59.1° and 
N c =6.1, yet again giving almost round figures. 

Assuming therefore that the flow is attempting to 
achieve a regular stable system with N c — 6 precisely, 
and taking q=2, we obtain from (36/35) Td z = 0.210 and 
Ps = 1.300, close to the observed but intermittent P3, but 
requiring a slightly larger r dz . Thus we may interpret the 
Q-mode as a failed attempt by the pulsar to switch to a 
second harmonic state with (from (26)) a particle drift of 
exactly half the pulsar rotation rate. It fails because such 
a drift is not consistent with 20 nodes, and requires them 
to be reduced. So why doesn't the pulsar simply reduce the 
number of nodes to a compatible figure? Possibly because 
the current flow from and to the polar cap, presumably pro- 
portional to the number of nodes, has to be maintained at 
some fixed rate whatever the level of the electric potential. 

Note that in the catalogue of 'drifting' pulsars of 
Rankin(1986) there are no fewer than 6 examples of pulsars 
with a P3 of 2.1-2.2 (the alias of p3 = 1.85), out of a total 
of 28. This seems a very high proportion, and it is possible 
that some - even all - these pulsars have exactly the same 
P3. In addition to PSR0943+10, the catalogued pulsars are 
PSR2303+30, PSR0834+06, PSR2021+50, PSR2310+42 
and PSR2020+28. PSR1933+16 has also been reported to 
show this phenomenon (Oster et al., 1977, Wolszscan, 1980)) 
and PSR1632+24 is a clear example (Hankins & Wolszczan, 
1987). It seems there is evidence of a pulsar subpopulation 
with this on-off property (see Fig 4 of Rankin 1986). Clearly 
aliasing bedevils a proper analysis of the drifting behaviour 
of many pulsars, and more careful investigation is needed 
and could provide a useful test of this model. 




:y[indcc 



Figure 5. The emission geometry for pulsars at a low angle of 
inclination. Note that although the rings of nodes both near the 
pusar surface (dark ring) and their mirrors on the null surface 
(finely dotted ring) encircle the magnetic axis, only the mirror 
ring at also includes the rotation axis. A low node drift (ie low 
subpulse drift-rate) therefore implies near corotation of the outer 
nodes on the surface of the corotating dead zone. 



7 THE INCLINED ROTATOR 

Many of the relatively stable (and slow-drifting) pulsars 
appear to have low anges of inclination (see the examples 
listed in Section 5.4 and the tables of Rankin 1993b). This 
is unlikely to be a coincidence. Near-axisymmetric geome- 
try means that quantities such as rj z , the distance between 
the outer nodes and inner nodes, and its associated time 
scale, t^ z , remain roughly constant along the trajectory of 
the nodes. Thus conditions for pair-creation will not vary 
too much as the nodes drift. Fig 5 shows graphically how 
that, although the observed nodes circulate the magnetic 
axis, the mirror nodes enclose both the magnetic and the 
rotation axis. The slow drift in pulsars of this type, often 
observed to be many orders of magnitude slower than the 
rotation period, therefore implies that the locations of the 
mirror nodes are almost precisely corotating with the dead 
zone surface on which they sit, although the real particle 
flows which produce them will be far from corotating. 

In near-axisymmetric conditions it is easy to see how 
drifting subpulses and the resulting emission cones can be 
maintained. But the model here, in contrast to the RS 
model, requires interaction between widely-separated re- 
gions of the magnetosphere, and it cannot be assumed that 
conditions for the creation of emission cones can persist for 
larger values of a. 

Figs 6, 7 and 8 show the fiducial plane which contains 
both the rotation and magnetic axes for increasing angles 
of inclination. For simplicity, at all angles a pure dipole ge- 
ometry is assumed rather than the more exact, but more 
complex, Deutsch (1956) solution for a rotating dipole in 
vacuo (Arendt & Eilek, 2001). Indeed an even better approx- 
imation may be the complex fully-corotating field derived 
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Figure 6. Connectivity at a = 35.3°. At this inclination the 
'jet' meets the light cylinder at the same colatitude as the light 
cylinder mirror point, enabling interaction between these regions. 
Furthermore the dead zone osculating point is capable of interac- 
tion with the opposite light cylinder mirror point. 




Figure 7. Connectivity at a = 54.7°. At this inclination the 'jet' 
colatitude coincides with that of the osculating point of the dead 
zone. For a larger than this the possibilities of regional interaction 
become greatly reduced 

analytically for inclined geometries by Beskin et al (1993). 
Nonetheless our assumption has the virtue that the null sur- 
face still intersects the plane in two straight lines (as in the 
aligned case), but now makes differing angles with the dipole 
axis. Here we focus attention on the changing locations of 
the mirror points in this plane as a increases, and their pos- 
sible interactions with other features in the azimuthal plane 
(the plane perpendicular to the rotation axis) . 

One such feature is a polar 'jet': It is known that most 



Figure 8. Connectivity at a ~ 90°. At near-perpendicular ge- 
ometries interaction is possible between opposite poles through 
the opposing deadzone osculating regions, and even through jct- 
to-jet contact 

pulsars display activity in the core region of their profile, 
so we can reasonably assume that a stream of particles em- 
anates from or infalls onto the magnetic pole along field-lines 
in a close bundle around the polar field-line. If the pulsar is 
not aligned with its rotational axis, this stream will cross or 
'impact' the light- cylinder at an angle which for simplicity 
we will assume to be a. However, near the light-cylinder the 
spinning of the star will twist the polar field-lines backwards 
in azimuth, together with the streaming particles, forming 
a ring-like structure at this colatitude. 

A second possibility is the osculating boundary of the 
dead zone at the light cylinder, represented in Figs 6 and 
7 by the extremes of the dotted lines. As is pointed out in 
Mestel (1999), particles constrained to rigidly corotate close 
to the light-cylinder can be expected to achieve sufficient 
inertia to drift across field-lines. These particles would rotate 
more slowly than the rigid rotation and form a ring around 
the rotation axis at the appropriate colatitude, encountering 
open fieldlines and possibly contributing to screening. 

As a increases from zero, the angle between the rota- 
tion axis and r* Cil , decreases monotonically, until a point is 
reached where this angle is itself equal to a. This is at 

a = arctan (4=) = 35 - 3 ° ( 43 ) 
v2 

At this angle the 'impact point' of the core jet is in the 
same azimuthal plane as one of the mirror points of the 
inner cone (see Fig 6). In the same configuration it can be 
seen that other mirror point of the inner cone (at r* c j) is at 
almost the same angle (= 2a = 70.6°) to the rotation axis as 
the osculating point of the dead zone(= 66.8°) and provides 
a second opportunity for azimuthal interference, this time 
between the second mirror point in this plane and, via the 
osculating point, the inner conal structure of both poles. 
At an inclination angle of 

a = arctan V2 = 54.7° (44) 
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a further geometric coincidence occurs (Fig 7) where the 
magnetic axis and the dead zone contact are at the same co- 
latitude on opposite sides of the rotation axis. This suggests 
a direct interactive link between the polar jet and the outer 
cone near the light-cylinder in the azimuthal ring perpendic- 
ular to the rotation axis. The locus of the inner cone mirror 
points at the light cylinder must also cross this ring at two 
places, so both cones and the core are capable of interacting 
in a single system. This may account for the complex be- 
haviour of PSR 1237+25 (Hankins & Wright 1980), a pulsar 
which is thought to have an inclination angle of 53° (Rankin 
1993b), and in which the behaviour of the inner and outer 
cone are clearly coordinated. 

The above picture of interactivity, based only on a ge- 
ometric analysis of the relative locations of the null surface 
intersections and the magnetic axis, lends a surprising de- 
gree of support to the results of RI, in which Rankin as- 
sembles two histograms of the numbers of each profile type 
versus a in 10-degree bins. As pointed out in Section 2, 
these results require extraneous - but not unreasonable - as- 
sumptions concerning the apparent angular size of the core 
region, yet they display a dependency on a which strongly 
supports the interactive picture developed here. The first 
(Fig 2 of RI) shows the histogram for pulsars with a promi- 
nent core component (type St), a young population with a 
mean age around 10 6 yrs. This reveals a sharp peak at 35°, 
which tails off up to 65°. In a further comment (in RII) 
on the distribution of type St pulsars, Rankin finds that St 
pulsars can be divided into two sub-populations, one having 
no cones at any frequency and the other developing inner 
cones at higher frequencies. The only factor discriminating 
between these populations is that the former cluster around 
35° and the latter around 50°. 

Using our interactive model we can interpret these ob- 
servations to mean that in young pulsars with narrow light- 
cylinders the core component only appears when it is pos- 
sible for the core jet to azimuthally interact, via a non- 
equatorial ring structure, with the mirror regions of the in- 
ner cone, and that that it is further enhanced by the in- 
ner cone's link to the dead zone (and hence to the outer 
cone), a = 35° therefore represents a 'resonant' inclination 
which encourages a strong jet, yet inhibits the formation 
of cones. For larger a, azimuthal interaction between the jet 
and the mirror regions remains possible in planes other than 
the fiducial plane, but will lack the 180° symmetry as the 
pulsar spins. As we approach the inclination where the jet 
and the dead zone can interact directly (a — 54.7°) an inner 
emission cone is able to form, suggesting that the function of 
its mirror region in the particle flow has changed now that 
the jet can interact with the opposite pole. 

The second histogram (Fig 4 of RI) relates to the 
broader class of pulsars with multiple cones or single outer 
cones (M and T types) and also shows a strong peak at 
35°, but with a somewhat wider spread. This suggests that 
interaction between different regions of the magnetosphere 
encourages emission cones to form for all pulsar types. Az- 
imuthal magnetospheric interactions with a ring structure 
enabling the current to close may indeed be a prerequisite 
for the 'pulsar phenomenon' to appear at all. But in these 
generally slower pulsars with wide light cylinders and weak 
jets (if any), the phenomenon takes a different form. Al- 
though the 35° peak strongly suggests interaction with the 



jet region is critical, the weakness of the jet no longer in- 
hibits cone formation. Cone formation now depends on a 
combination of the altitude of the mirror regions operating 
in a non-axisymmetric geometry, the azimuthal position of 
the jet in relation to the mirror regions, and almost cer- 
tainly the level of activity in the jet. The jet and the cones 
will exist in a symbiotic relationship. 

For a between 60° and 85° the number of pulsars of 
all conal types falls off dramatically in the histograms. At 
these angles the mirror points on the light-cylinder, the dead 
zone contact and the jet are now widely separated from each 
other in colatitude, discouraging interaction (indeed any in- 
teraction would probably be between the opposite poles). 
Furthermore the field lines linking the surface to the mirror 
points for the outer cones (on the surfaces of the dead zone) 
are becoming severely distorted, demanding a long and com- 
plex route across the rotation axis on the one side, and with 
the distance to the mirror point (r^ z 2 = 0.06i?; c ) on the 
other side diminishing well below the levels demanded by 
the criterion of the previous section. Similar distortions can 
be seen in the inner cone mirror points, and recall the result 
of the previous section that, for whatever reason, extreme 
high or low mirror points appear to inhibit cone formation. 

But finally, at inclinations close to 90° (see Fig 8), the 
situation changes dramatically and a strong peak appears 
near a = 90° in both histograms (see again Figs 2 and 4 of 
RI). Now the dead zone makes contact on both sides of the 
light cylinder and fills an azimuthal ring about the rotation 
axis. Interaction, now between opposite poles, is again es- 
tablished and pulsars of all types are indeed found to cluster 
at this angle. Direct contact between the polar jets is quite 
possible and and is supported by examples of observed in- 
terpole interaction in a number of pulsars (Fowler & Wright 
1982, Gil et al 1994, Biggs 1990). 

Although we have compromised with reality by con- 
sidering only a strictly dipole geometry, the results of this 
section show that it is possible to account for the observed 
dependence of pulsar type on a by invoking light-cylinder in- 
teractions as the catalyst for pulsar emission. It would seem 
plausible that the observed distribution of pulsar inclina- 
tions is not an evolutionary effect, and that therefore many 
rotating magnetic neutron stars, especially at high inclina- 
tion, fail to become observable pulsars because they lack the 
necessary connectivity. 



8 THEORETICAL ISSUES 

We have now arrived at the point where we must consider 
what physical system might lie behind the essentially geo- 
metric picture we have built so far. 

8.1 The global magnetosphere 

The mirror points we have identified are locations where 
transitional behaviour in the charge flows of the magneto- 
sphere might be expected: At the light-cylinder intersection 
the particles will have to negotiate the sign-change in the 
charge density, and also somehow keep their velocities sub- 
luminal as they cross the light- cylinder. At the dead zone in- 
tersection, there must be a sharp lateral transition layer be- 
tween the corotating region and a subrotating flow of charge 
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moving along the field-lines into and out of the gap. It is 
therefore not unreasonable to suspect that these two regions 
may influence, or even be responsible for, the cap emission. 

The significance of these intersections was anticipated 
in the axisymmetric magnetosphere model of Mestel et al. 
(1985, henceforth MRWW), and similar ideas are present 
in that of Shibata (1995), designed to explain gamma-ray 
producton in fast pulsars. Both papers are based on the 
principle that to create the torque necessary to brake the 
star, the current must close outside the light-cylinder, leav- 
ing the magnetosphere at a higher latitude and returning at 
a lower. MRWW envisage electrons streaming from the po- 
lar cap to the light-cylinder intersection, at which point they 
are forced to leave the field-lines, and postulate the creation 
of a gap void between this point and the dead zone. Thus 
the return current must find its way back between the outer 
gap and the corotating dead zone. There are echoes of this in 
the configuration described earlier for a = 54.7° (see Fig 7), 
although the axisymmetric nature of the theoretical models 
inevitably precludes the possibility of using azimuthal flow 
to complete the current circuit. Such connectability within 
inclined rotators may turn out to be a key ingredient. 

8.2 The feedback process 

In requiring a two-way flow between the polar cap and the 
gap region on both the critical fieldlines, the model here also 
suggests a feedback process. The simplest model would ac- 
celerate particles from the cap all the way to the null line 
without completely screening the electric field parallel to 
the magnetic field lines, and would there create pairs in an 
outer gap. Some positrons would then be available to com- 
plete the feedback loop implied by the 'mirror' concept and 
return to the surface. Somehow the returning positrons must 
interact with the electron flow, and so on. Models relevant 
to our purpose are those which are set in conditions of near 
GJ screening (so-called space charge limited flow models) 
and include the accelerator models of Arons et Scharlemann 
(1979), Mestel & Shibata (1994), Shibata (1997) and Mes- 
tel (1999). The attraction of these models here is that the 
outflowing particles (envisioned by these authors as elec- 
trons) are only gradually accelerated to high altitudes be- 
fore Lorentz factors of sufficient magnitude to initiate pair 
production are achieved. 

Recent promising developments of the space-charge- 
limited flow models (Arons & Scharlemann 1979, HAa,b, 
HM, HMZ), recognising the the role of (nonresonant) in- 
verse Compton scattering of thermal photons - rather than 
the curvature radiation operating in young pulsars (HAa, 
Harding & Muslimov 2001) - as the trigger for pair produc- 
tion in 'slow' pulsars, demonstrate that the created pairs 
cannot fully screen the electric field parallel to the mag- 
netic field. This feature (described in HM) fits well with the 
requirement here that particle acceleration should continue 
up to the outergap location, and is further supported by 
the critique by Jessner et al (2001) of dense pair production 
models. However space charge- limited models rely heavily on 
suitable axisymmetric dipole geometry to drive the acceler- 
ation (by means of divergencies from GJ charge densities). 
Observations of drifting subpulses, and the explanation for 
these offered here, imply a highly non-axisymmetric poten- 
tial and charge distribution above the pole (although still 



with a dipolar magnetic field) . It would be of interest to ex- 
amine the consequences for these models of, say, a regular 
azimuthal variation in potential. 

In order to create a feedback to the surface we also 
need a model for pair creation in the outer magnetosphere 
to follow the magnetosphere acceleration, and one which will 
work in 'slow' pulsars, i.e. those with periods of the order 
of a second. The model of Shibata (1995) - designed for fast 
pulsars - stressed that any supposed inner accelerator sys- 
tem must influence the electrodynamics of the light-cylinder 
region in order to successfully dispose of the star's angular 
momentum. The analysis here suggests that the accelerator 
is fully integrated into the dynamics of the outer magne- 
tosphere, even for pulsars which are not observed to emit 
gamma-rays, and that pulsars have found a way of feeding 
back changes in the outer magnetosphere to the polar cap. 
As yet, the possibility of weak, non-cascading, intermittent, 
and also non-axisymmetric, pair production in an outergap 
acceleration zone has not been examined, and it therefore re- 
mains far from clear whether the outergap can sustain itself 
in its highly diminished role. But the detailed work of HSa,b 
(also Hirotani 2000), although focussed on high-energy pair 
production, intriguingly suggests that the location of the 
outergap may depend on the inflowing current from the star 
and hence be part of the feedback process - thereby hinting 
at a mechanism for mode-changing. 

8.3 The emission process 

In the pair-creation cascade at the mirror points the elec- 
trons will be accelerated and injected into the wind zone 
beyond the light-cylinder. But, as Shibata (1997) noted, at 
least a small fraction of the positrons will trail behind the 
cascade process and feel the same potential as the primary 
electrons and return to the surface. As they approach the 
star ( e.g. Michel 1992), the positrons will acquire a 7 factor 
sufficient to pair-create and an avalanche of bunched pairs 
will be formed some 10 star radii above the surface. This 
process is analogous to that of cosmic ray entry into the 
atmosphere, as noted by Melrose (1996). In Michel's model 
the downward coherent radiation produced by this process 
is reflected from the surface. The radiation is exactly ra- 
dial, has the dimensions of the polar cap (and reverses its 
polarisation), and would seem to be well-suited to explain 
the emission features and dimensions of the pulsar profile's 
central core. Rankin (in RI) observed that the central (core) 
component of pulsar profiles scales precisely with the open- 
ing angle of a dipole polar cap at the surface of the star, 
just what would be expected from Michel's model. Effec- 
tively, the reflected radiation will be a miniature image of the 
conal structure, and the time lag behind the conal radiation 
at 220km will shift the core position by around 0.7ms. It is 
important to realise that the number of returning positrons 
must be well below the GJ value so as not to 'poison' the 
outflow (Michel 1991, Lyubarskii 1992), and will impact only 
a fraction of the polar cap surface. Otherwise, as Shibata et 
al (1998) have pointed out, the polar cap would emit X-rays 
above the observed limit formulated by Becker & Truem- 
per(1997). 

On the other hand, X-rays from the surface may be part 
of the mechanism which produces pairs in the outer gap 
(Harding 2000, Romani & Yadigaroglu 1995). This might 
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explain why a cone cannot be formed when linked to a gap 
above a fixed altitude, rather than to a fixed proportion 
of the light-cylinder radius. Similarly, below a certain fixed 
altitude pair-production, stimulated by surface X-rays, may 
result in so many backflowing positrons that the flow is 'poi- 
soned' and becomes unstable (Lyubarskii 1992). 

The pairs avalanche initiated by the backflowing 
positrons will clearly heat the surface immediately below 
it and release electrons by a thermionic process. These new 
primary electrons will have a low 7 and must pass through 
the oncoming avalanche above it. If they are to emit the 
coherent radiation seen in pulsar profile cones then it must 
be assumed that they are bunched through instabilities in 
the strictly layered avalanche. This configuration does not 
seem to have been analysed in the literature, but Melrose 
(1996) has pointed out that such a model avoids many of the 
difficulties (nothing to maintain a pancake structure, emis- 
sion angles too narrow (-)) associated with models based 
on purely outflowing streams. 



9 CONCLUSIONS 

In this paper an attempt has been made to create a feedback 
model for the pulsar magnetosphere which can be directly 
related to observations. It sees the magnetosphere as an inte- 
gral whole, where the polar cap and the outergap are mirrors 
of each other. It is suggested that the inner and outer cones 
found in integrated profiles reflect the two intersections of 
the Holloway null surface with the light cylinder and the 
corotating dead zones respectively, and the core components 
of profiles quite literally reflect the feedback from these in- 
tersections. This simple interpretation of the profiles can 
reproduce with some accuracy the observed dependence of 
profile features on period and angle of inclination. The ob- 
served regularity of profile features from pulsar to pulsar is 
difficult to explain by the arbitrariness of complex surface 
magnetic fields. 

An important feature of the model is that for the first 
time it suggests an natural evolutionary link between old 
slower radio pulsars and young Crab-like pulsars whose pre- 
dominantly high-energy X-ray/7-ray emission is thought to 
emanate from the outergap regions of an inclined rotating 
dipole (Romani & Yadigaroglu 1995). As a pulsar ages, cas- 
cades in the outergap accelerated regions fade (Ruderman 
& Cheng 1988), but a meagre yet critical pair production 
still lingers, which, by interacting with polar cap produc- 
tion, creates the coherent radio emission we observe in the 
complex radio profiles. 

The drift of particles in the entire magnetosphere often 
gives rise to discrete rings of emission nodes near the po- 
lar cap. They can move in either direction, and will have 
differing behaviour in different rings. The particle flow and 
the observed motion of the nodes may be chaotic, or adopt 
a 'resonant' state where the particle circulation is harmon- 
ically coupled to the star's rotation rate. PSR0943+10 is 
shown to be in such a state, and it is argued that conditions 
in other stars may permit the star to cycle through several 
resonant states of stepwise increasing charge screening along 
open field-lines. In general, the observed subpulse drift pat- 
terns are only the residual drift of the particle circulation, 



and hence are sensitive to small variations in the ambient 
electrical potential. 

The angle of inclination plays an important role in 
determining whether a pulsar successfully creates emission 
nodes and hence observable coherent radiation. It changes 
the position of the null surfaces and thereby fixes the de- 
gree to which the magnetosphere is interconnected by means 
of rings of azimuthally drifting particles around the light- 
cylinder. The nature of the observed dependency of incli- 
nation on conal type suggests that the flow of particles be- 
tween the core region of the magnetic polar cap and the 
light-cylinder is an integral part of the system, and this can 
be shown to further imply that the two magnetic poles will 
generally be linked in a single system. 

A further constraint on the formation of cones, sug- 
gested by the observed dependency of conal type on the 
rotational period and made plausible by the interpretation 
of this model, is that emission nodes, and their resulting in- 
tegrated emission cones, can only form if the mirror points 
lie between fixed altitudes of approximately 20,000km and 
70,000km. 

If core emission is caused by radiation reflected from 
the neutron star's surface, then it is a prediction of the 
model that on close inspection the core component will have 
a structure of miniature cones (an effect possibly easier to 
detect in younger pulsars). 

However the model leaves many theoretical questions 
unanswered. How exactly can the pair creation process be 
made to work in the outer gap (HSa,b), and how many 
positrons can really be available for the backflow to the sur- 
face? To what degree is the proposed slow acceleration of 
particles between the polar cap and the outergap (HAa,b, 
HM, HMZ) affected by a returning charge flow of opposite 
sign and with azimuthal dependence? Can an emission pro- 
cess be made to work as slow electrons pass upwards through 
the descending pair avalanche? A more radical question is 
whether we need particles from the surface at all to make the 
radio emission. Then the model could apply to pulsar and 
antipulsar alike. Finally we could speculate, as in the global 
model of MRWW, that the outer ring represents a net flow 
of negative charge to the surface, surrounding and balancing 
the positive flow within the inner cone, and thereby solving 
the long-standing current balance problem, although in suf- 
ficiently inclined pulsars non-equatorial rings could close the 
current in the azimuthal plane (see Section 7). This model 
has many features needing theoretical attention, but its aim 
is to forge a closer link between observation and theory, and 
the author will be pleased if something of this is achieved. 
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